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Abstract. We show that one can define through the symmetry approach a procedure to check the 
hnearizabihty of a difference equation via a point or a discrete Cole-Hopf transformation. If the equation 
is linearizable the symmetry provides the linearizing transformation. At the end we present few examples 
of applications for equations defined on four lattice points. 



1. Introduction 

In recent years the use of partial difference equations has been playing an increasing role in physics and 
in mathematics. From one side discrete systems are believed to be at the base of the basic laws of physics 
(see for example the recent literature on quantum gravity [22]) and on the other side with the increasing use 
of computers, discretizations are playing an increasing role for solving numerically differential equations [7]. 

Calogero [S] introduced a heuristic distinction between nonlinear partial differential equations which 
are " C-integrable" and " S-integrable" , namely, equations that are linearizable by an appropriate change 
of variables (i.e. by an explicit redefinition of the dependent variable and maybe, in some cases also 
the independent variables), and those equations that are integrable via the Inverse Scattering Transform 
(1ST) 4 . In the cited reference Calogero used the asymptotic behavior to find equations belonging to these 
two classes of equations as it usually preserves integrability and linearizability. However this approach, 
even if very fruitful is often very cumbersome and not always exhaustive (see for example the results of [S] 
for discrete equations). 

A more intrinsic approach is based on the existence of symmetries. The well-known notion of higher 
symmetry is at the base of this approach together with the notion of 1ST. The mutual influence of these 
theories has led to the fundamental abstract concept of formal symmetry PSII^ , more basic than symmetry 
as it provides also higher conservation law, Backlund transformation and Lax pair representation. In this 
sense a formal symmetry is a universal object. However the derivation of formal symmetries is difficult 
to apply in the case of difference equations _2Cl, . Moreover the distinction between C and S integrable 
equations is only at the level of the conservations laws as linearizable equations have no local conservations 
laws of arbitrary high order [55]. 

A well established result in the framework of Lie theory for proving the linearizability of nonlinear 
Partial Differential Equations (PDEs) is provided by Kumei and Bluman [T^ (for a recent extended review 
see 1]) based on the analysis of the symmetry properties of linear PDEs. Following the analogy of the 
continuous case we will formulate a similar theorem for linearizable Partial Difference Equations (PAEs) 
by which we recover part of the results obtained by assuming the existence of linearizing transformations 
[TH[Tni[25- Partial results in this direction for the case of difference equations defined on a fixed lattice 
have been obtained by Quispel and collaborators [51 1^1^ . 

In Section [5] we review for the sake of completeness and clarity the results presented by Bluman et. 
al. in [T] for the linearization of nonlinear PDEs both by point and nonlocal transformations and consider 
examples of their applications. Section[3|is devoted to the introduction of PAEs, their symmetry analysis 
and the discretization of the theorems presented in [1 . A series of examples is presented in Section [5] 
while in Section [S] we summarize the results and present conclusions and outlooks. In an Appendix a 
theorem complements the results of Section [3] by showing that equations possessing an infinite dimensional 
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symmetry algebra depending on an arbitrary function are linear and the arbitrary function satisfies a linear 
homogeneous equation. 



2. Linearization of PDEs through symmetries 

In |13j Bluman and Kumei introduce a series of theorems dealing with the conditions for a nonlinear 
PDE to be transformable into a linear one by contact transformations. Here, in the following, we will 
limit ourselves to the case of just point transformations as these are the relevant ones in the discrete 
case [TTKIH]. In more recent works the same authors [5] extended the consideration to the case when we 
have non-invertible transformations between a nonlinear and a linear PDE. 

The basic observation is that a linear PDE, 

Cv{y) = F{y), (1) 

where £ is a w-independent but possibly y-dependent linear operator and ^{y) is the inhomogeneous term, 
has one point symmetry of infinitesimal symmetry generator 

^ = (2) 

depending on a function w which satisfy the homogeneous equation 

Cw{y) = 0, (3) 

as any solution of ([T]) is the sum of a particular solution plus the general solution of the associated 
homogeneous equation. As the existence of an infinitesimal generator of the form ^ is preserved when 
we transform a linear equation into a nonlinear one by an invertible point transformation, following [T] 
Section 2.4, we can state in the following the conditions for the existence of an invertible linearization 
mapping of a nonlinear PDE: 

Theorem 1. A nonlinear PDE 

£n{x,U,Uj;,- ■ -Unx) = (4) 

of order n for a scalar function u of an r -dimensional (r > 2) vector x will he linearizable by a point 
transformation 

w{y)^f{x,u), y = g{x,u), (5) 
to a linear equation (0) for w if it possesses a symmetry generator 

r 

X = '^Ci{x,u)dx, + (t){x,u)du, ^t{x,u) = at{x,u)w{y), (6) 
1=1 

(t){x, u) = uj{x, u)w{y), 

with uj and Ui given functions of their arguments and w{y) an arbitrary solution of 

Following |13) we can state the sufficient conditions for the existence of an invertible linearization 
mapping of a nonlinear PDE. This theorem defines the transformation ([SJ: 

Theorem 2. // a symmetry generator for the nonlinear PDE as specified in Theorem{l\ exists, the 
invertible transformation which transforms ^ to the linear PDE (0j is given by 

yi = i = I,--. ,r, (7) 

w = '^'{x,u). (8) 

where ^i{x, u) are r functionally independent solutions, i — 1, ■ ■ ■ ,r, of the linear homogeneous first order 
PDE for a scalar function u) 

r 

ai{x, u)^{x, u)x^ + uj{x, u)^{x, u)u = 0, (9) 
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and '9{x,u) by a particular solution of the linear inhomogeneous first order PDE for a scalar function 

r 

at{x, u)'^{x, u)x, + uj{x, u)'if{x, u)u = 1- (10) 

i=l 

If a given linearizable nonlinear PDE does not have local symmetries of the form (jS]), i.e. its local 
symmetries do not satisfy the criteria of Theorem [TJ it could still happen, as shown in [T] Section 4.3, 
that a nonlocally related system has an infinite set of local symmetries that yields an invertible mapping 
of the nonlocally related system to some linear system of PDEs. Consequently, the invertible mapping 
of the nonlocally related system to a linear system will provide a nonlocal (non-invertible) mapping of 
the given nonlinear PDE to a linear PDE. This non invertible transformation will be a kind of Cole-Hopf 
transformation [HIII^]- In this case, however, we have to generalize Theorem [5] to take into account the 
fact that we are dealing with a system of equations. 

Theorem 3. Let us consider a system of nonlinear PDE 

£^^\x,U,V,U^,Va:, ■ ■■Unx,Vnx) = 0, £^^\x , U, V , U^, , ■ ■■Unx.Vnx) = (11) 

of order n for two scalar functions u and v of an r -dimensional (r > 2) vector x which possesses a 
symmetry generator 

r 2 

X = ^^^\x,u,v)dx, + 4'{x,u,v)du + ip{x,u,v)dv, Cix,u,v) ^^^aj{x,u,v)w''^'> (y), (12) 
i=i j=i 

2 2 

(l){x,u,v) = ^I3j{x,u,v)w^^\y), i(;{x,u,v) =^'yj{x,u,v)w^^\y), 

with aj, (3j and "fj given functions of their arguments and the function w — {w'^^\y) , w'^'^^y)) satisfying 
the linear homogeneous equations 

M{,y)w{y)^Q, (13) 

with y an r -dimensional vector depending on u, v and the vector x and Ai is a 2x2 matrix linear operator. 
The invertible transformation 

w(i)(y) = i^(i)(a;,M,w), ^ f(2)(x,u,w), y = G{x,u,v), (14) 

which transforms Ul\) to the system of linear PDEs il3\) is given by r functionally independent solutions 
Gi{x,u,v) with i = 1, • • • ,r of the linear homogeneous first order system of PDEs for a scalar function 
G{x, u, v) 

r 

al{x, u, v)Gx, + (ik{x, u, v)Gu + ^k{x, u, v)Gy = (15) 

i=l 

and by a particular solution of the linear inhomogeneous first order system of PDEs for the function 
F = {F^^'>{x,u,v),F^^'>{x,u,v)) 

r 

^k{x, u, v)F^ + I3k{x, u, v)Fi^^ + 7fe(a;, u, v)F^^^ = (16) 

i=l 

with (5^ the standard Kronecker symbol. 

For the sake of completeness and to clarify the application of the theorems presented above, in view 
of the discretization which will be presented in the following section, we consider here one example of 
linearizable nonlinear PDEs belonging to each of the two cases presented above. 
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2.1. A nonlinear PDE linearizable by a point transformation. It is well know, see for example 
Olver book [21], that the potential Burgers equation 

Ut = U^x + {UxT, (17) 

is linearizable by a point transformation. In fact the infinite dimensional part of the infinitesimal generator 
of its point symmetries is given by 

X ^v{x,t)e-'^du, (18) 

where v{x,t) satisfies the homogeneous linear heat equation vt — v^x — 0. 

The conditions of Theorem [T] are satisfied with lo = e~" and ai = 0. We can apply Theorem [3] and we 
get = X and $2 = ^ as from ([1]) = while from ([TU)) '^{x, u) satisfies the equation '^{x, u)u = e" i.e. 

u = \og^{w). (19) 

Eq. ([T5)) is the linearizing transformation for the potential Burgers equation ([TT]) . 

2.2. A nonlinear PDE linearizable by a non invertible transformation. The standard example in 
this class is the Burgers equation 

Ut^Uxx-UUx^[Ux-]^u'^]x, (20) 

linearizable by a Cole-Hopf transformation. As (j20p has no infinite dimensional symmetry algebra but it 
is written as a conservation law we can introduce a potential function v{x,t) and (|20p can be written as 
the system 

vx = 2u, (21) 
Vt = 2ux - y?- 

Applying Theorem [3] we can find an infinite dimensional symmetry of the form of (jl2p . In fact, solving 
the determining equations, apart from terms corresponding to a finite dimensional algebra, we obtain an 
infinite dimensional dilation symmetry given by 

tp — 4:Uj{x,t)e^ , (f) = -ipx + ut/jy — e~[2uJx + uju] (22) 

where uj{x, t) satisfies the linear heat equation uit — ^xx = 0. 

The linearizing transformation can be obtained from Theorem [3] Let us define w^^^y) = Lij{x,t), 
w^'^\y) = uJx[x,t) and take as functionally independent solutions of (fT5)) Gi — x and G2 = t. As = 
and 71 = 4e^, 72 ue^ and /32 — 2e^, we get as a particular solution of (fTB]) 



F^^^^-e^^, i^(2) = iue~^. (23) 

Eq. (|14l) implies uj — — e^^ and uJx — ■^'ue~^ and from it we obtain as a linearizing transformation the 
Cole-Hopf transformation 

^^x 

u — —2 — . 

w 

3. PAES AND THEIR LINEARIZATION 

We will consider here the problem of the linearization of a partial difi'erence equation. 

As is well known difference equations can be obtained on a given grid, starting from physical or chemical 
or biological problems dealing with lattice systems, or as discretization of differential equations with 
symmetries which we want to preserve. In this second case we can use the freedom of the grid to get a 
difference scheme which preserve part of the symmetries of the continuous equation. 

In this first approach to the problem of linearization of difference equation we will consider the case 
when the grid is preassigned and assumed to be fixed, with constant lattice spacing. Moreover for simplicity 
we will consider autonomous equations defined on a two dimensional grid so that there is no privileged 
position and we can write the dependent variables just in terms of the shifts with respect to the reference 
point Un,rn — uq^q ou the lattice. 
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A partial difference equation (PAE) of order N ■ N' for a function Un,m will be a relation between N ■ N' 
points in the two dimensional grid, i.e. 

£n-n' (uq^o, ui^o, ■ ■■ , UN,a, uq^i, • • • , un,i, • • • , un,n'^ = 0. (24) 

A continuous symmetry for equations of the form (j24l) . where the lattice is fixed, i.e. the two independent 
variables Xn,m and tn,m are completely specified as x„ = h^n + Xq and tn_m = htm + to with h^, ht, xq 
and to given constants, is given just by dilations 

It is easy to show that a linear PAE of orderA^ • A^' for a function Vn,m 

{N,N') 

Tn-n' =b{n,m) + ^ a'■'^n,m)vn+i,m+j = 0, (26) 
(^j)=(o,o) 

has always the symmetry 

where (j)n,m is a solution of the homogeneous part of ([25]) . It is not at all obvious, however, that an 
equation (|24p having a symmetry (1271) is linear and that the function (j)n,m must satisfy a homogeneous 
linear equation. We leave to the Appendix the proof of this proposition. The symmetry (j27p is due to 
the superposition principle for linear equations. If the nonlinear equation (|24p is linearizable by a point 
transformation then the symmetry (I27p must be preserved. This is the content of the Theorem [T] we 
presented in the previous Section in the case of PDEs and this must still be valid here. So we can state 
the following theorem: 

Proposition 1. An autonomous linear PAE \24^ is linearizable if it has a point symmetry of the form 

where the function (f>n.m satisfies a linear PAE of the form i26]) with b{n,m) = 0. 

The proof of Theorem [1] and of the following Theorem U] are the same as those presented by Bluman 
and Kumei [2] in the continuous case and so we will not repeat them here. 

As in the case of PDEs we can present the following theorem which provide the transformation which 
reduces the equation to a linear one. In this case, as the independent variables are not changed, the 
transformation is given just by a dilation. So we have: 



Theorem 4. The point transformation which linearizes the nonlinear PAE 
is obtained by solving the differential equation 

an,m{Un,m) = 1- (oO) 

As in the continuous case, if (|24l) has no symmetries of the form considered in Theorem [T] we can 
introduce some potential variables. On the lattice there are infinitely many ways to introduce a potential 
variable as there are infinitely many ways to define a first derivative. Thus it seems to be advisable to 
check the equation with a linearizability criterion like the algebraic entropy [26] before looking for potential 
variables. 

The simplest way to introduce a potential symmetry is by writing the difference equation (|24p as a 
system 

In such a way 



N-N 



However it is easy to show in full generality that the symmetries for (PT|) and for are the same. 
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We can introduce potential symmetries by the following system, 
In such a way 

Sn-N' (uofi, -Ui,o, • • • , UNfi, Uo,l, • • • , UN,1, • • • , UN,N') = ^ ~ [^+1,™ " £L%]- (34) 

i.e. the nonlinear difference equation is written as a discrete conservation law. As p3p is a system, to 
construct the symmetries we have to generalize the linearization theorem as we did in the continuous case. 

Theorem 5. Let us consider a system of nonlinear PAEs 

"^ri.m i^n^m-j ' ' ' ; ^n^va ; ' ' ' ) ^n rn (^n,m i ' * ' i ^n^va ; ' ' * ) — (35) 

of order N ■ N' for two scalar functions Un^m o,nd Vn.m of two indices n and m which possesses a symmetry 
generator 



n.ni 1 



with /3^^' and 7*--'-* given functions of their arguments and the function w — (wi^L,wi^L) satisfying the 
linear equations 

^n,m^n,m — 0- (37) 

The invertible transformation 

which transforms i35\} to the system of linear PDEs |g7| j is given by a particular solution of the linear 
inhomogeneous first order system of PDEs for the function K = {Kn.ln{un,rmVn,m 

/^n,m(^n,ni; '^n,m)K^^^ + ^n"Jni'^'ri,rm '^n,m)K^^^ — S'j^, (39) 

where S-^, is the standard Kronecker symbol. 

4. Examples. 

Here we present a few examples of linearizable PAEs. For concreteness and for comparing with the 
previous literature [15j we limit ourselves to the case when the nonlinear difference equation involves at 
most 4 lattice points. 

4.1. Classification of PAEs on a square lattice linearizable by point transformations. We con- 
sider a general autonomous equation defined on a square lattice: 

J"(mo,o, uo,i, ui,o, ui,i) = 0, (40) 

where, for convenience, we just write down the shift with respect to the reference point of indices n, m. If 
we assume that ui i is present in (|40p then we can rewrite the equation as 

= F{uq,q, uo,1i "i.c)- (41) 
Following Theorem [1] we look for an infinitesimal symmetry generator of the form 

Xqm ^ ao^oiuo,o)4>o,oduo,oi (42) 
where the function </> solves a linear homogeneous equation, i.e. 

C(l)o,o - 0, C = a + bTi + cT2 + dTiTs , (43) 
with Ti and T2 operators such that Ti(/!)o,o = 0i,o and 7200,0 = </'o,i- If d 7^ then we can write (|43l) as 



<^i,i = ~^[«<?^o,Q + &01.O + «/'o,i]- (44) 



In this setting uo.ii '^'j.Oj (^o.i and (f>jfl, with i,j = 0, 1 are independent variables. If (|42l) is a generator of 
the symmetries of (PT|) then we must have 

XJ" = o F„o 000, 0^0, 0(^^0,0) + -F',ui,o</'i,oai,o(wi,o) + -P',«o,i'?^o>i"o>i("o,i) = (45) 

= </'i,iai,i(wi.i) = -^[a'/'o.o + ^'/'i.o + c(/!'o,i]ai,i(-^("o,o, ""0,1' ""i.o))- 
a 

As 00,0; 01,0 a-nd 0o,i are independent variables, we obtain from (j45p three equations relating the function 
a, intrinsic of the symmetry, with the function F, intrinsic of the nonlinear equation: 

aai^i{F) + di^,„o.(,ao,o(wo,o) = 0, 6ai,i(F) + cfF„^ oai,o(ui,o) = 0, cai,i(F) + dF„o,iao,i(Mo,i) = 0. (46) 

As in (|46p . up to a constant, the first term is the same for all three equations, we can rewrite them as a 
system of PDE's for the function F depending on a 

--P',«o,o"o,o(mo,o) = -^F^ui_o(^i,o{uifi) = -F,uo,i(^o,i{uo,i), (47) 
which can be solved on the characteristic, giving as a function of the symmetry variable 

^ ^ ag{uQ,Q) + bg{ui,o) + cg{uQ,i), a{x) ^ — (48) 

Introducing this result in Theorem S] we get '0(^0, 0) — J^°'° gx{x)dx — g{uofl) + k, with k an arbitrary 
integration constant. Then (|46p gives that any linearizable nonlinear PAE on a four points lattice must 
be written as 

dF^^ + a{F{0) =0^ F = g-\^^), (49) 

where by g^^{x) we mean the inverse of the function g(x) given in (j48p . 

Let us notice that from (|47l) we can get the six linearizability necessary conditions we introduced in [T6] 
to classify linearizable, multilinear equations on the four lattice points, that is 

p 

A{x,ua^i) = — 

-^,«i,o 

D / \ • ."0,0 

B(x,ui^o) = 



= 0, 



a 

=111, o=a; ^ 


Vx, 


"0,1, 


(50a) 


a 

— U{) l—X ~; 

c 


Vx, 


"1,0, 


(50b) 


c 

= "0,1=2; — ^5 


Vx, 


"0,0, 


(50c) 


Vuo,o, "1,0, 


"0,1, 




(50d) 


Vuo,o, "1,0, 


"0,1, 




(50e) 


Vuo,o, ■"1,0, 


W0,1- 




(50f) 



C {x,Uq^o) = 

^,"1,0 

d F,uo^o _ 
duo,i F^ui^o 
d -P'.mq.o 
9"i,o ^","0,1 
d f,"o,i 

So linearizable equations on four lattice points are characterized by a function g{x) and its inverse. As 
a trivial example we can choose 5(2;) = and we get that the nonlinear equation ui.i = log(ae"°''' + 
^g"i,o _|_ryg"o,i _j_ linearizes to ipi^i = a^o,o + /3V'i,o + 7''/'o,i- The corresponding function a is a{x) = 
and the linearizing transformation is V'0,0 = e""-" + n. 

In |15] we have shown that there is a multilinear equation on the square lattice belonging to the Q"*" 
class which is linearizable. It is interesting to find the corresponding function g(x) in term of which 
we can linearize it. The function F in this case is a fraction of a second order polynomial over a third 
order polynomial. The only function g which provides this structure is g{x) — ^^^"'^^^ which gives F = 



C-?o 



where f = -„ — ^^-ra — h 7, n? — h 1 — ^-nr- In this case the linearizing transformation is 

^ £lU0.0+«0 «lMl,0+';o tlMO,l+*0 ° 



the linear fractional function ^"0,0 ("0,0) = £1 no ^0+^*0 ^ ^^^^^ arbitrary constant. 



_____ /i2\ 

4.2. Linearizable potential equations. For the sake of simplicity we set in £n,m = Uri,m- If we 

want the equation (p4|) to be on the square we have to choose Snjn — 9n,m{un.m, Un+i,m)- The application 
of the prolongation of the infinitesimal generator (pS)) to the second equation in gives [12] 

'/>o,o(mo,o,wo,o) = -00,1(^0,1, Wo, i) - ■0o,o(uo,o, wo,o), -> V'o,o(wo,o, vo,o) = V'o,o(wo,o)- (51) 
Then the prolongation of the infinitesimal generator applied to the first equation in gives 

V'i,o(wi,o) - ?Ao,o(wo,o) = [-00,1(^0.1) - ■0o,o(wo.o)]|^ + [V'i,i(«i,i) - V'i,o(wi,o)]f^, (52) 

where 

1^1,0 = Wo, + ffo,o(Mo,o, ""1,0), '^0,1 = t^o,o + "0,0, I'l,! = wo,o + "1,0 + .go,o(Mo,o, "i,o)- (53) 

To comply with Theorem [5] we look for an infinitesimal coefficient of the infinitesimal generator p6p of the 
form 

00,0 = w'o!d7o!o(wo,o) + w^o.hoM^ofi), (54) 
where the functions Wn ln and Wn m satisfy a linear partial difference equation on the square 

'^0,0 = ao!o"'M +4^ow[^o +4^o^t'i^i, (^^) 

y,(2) _ ^(1) (2) ,(2) (2) ^(3) (2) 

"'O.O ^ ^^0,0^0,1 + "o,o"'i,o + "o,o'"i,i- 
Introducing (|54l [55)) into (1521) and taking into account that we can always choose Wq^I, 'w[^l, w^^l, wl^l, 
w^^Q and w^j^l as independent variables we get the following system of coupled equations for 7^^^ 

7l!d(t'o,o+go,o)(l + fe)-4'd7^^,](t'o,o)(l + |f^) =0, (56) 

7^^^(^'o,o + "o,o)i^ + a^oho'Mvo,o){l + ^) =0, (57) 

l[]livo,o + "1,0 + 9o..o)t^^ + 4'd7d^,](«o,o)(l + =0, (58) 

and similar ones for the function jn}n{vQ,Q). Adding (1561) multiplied by a'j^Q to (j58p multiplied by Oq^q we 
get 

(3) (1)^ I I ^5o,o\ , (2) (1), , , ^9go,o „ 
^0 o7i ("0,0 + 50,0) (1 + ^ j + % o7i 1 (wo,o + "1,0 + 50.0) = 0, (59) 

an equation similar to (j57p . i.e upshifting by one the first index in (|57p and comparing the result with ([5 
we get a discrete equation for gn,m 



(3) /1 1 9go^\ 9go,o 

"•0,0 V am.o/ _ „(2) am.o 

(1) (i , 99MI^ " 

"1,0 V-^ ^ awi.o/ 9«i,o 



In (j60p appears the function gi^Q — (7i.o(ui,o, "2.0) a-nd if ^^^^'^ ^ ^ we get a linear differential equation 
for go whose solution is 50,0 — 5o^o("o,o) + 5o^o(wo,o)"i,o- Introducing this solution in (1501) we get go = 
So^o + 3o^d("o.o)"i.o +5o^d"o,o, i-e. a linear equation. By choosing, in place of (j55l) the most general linear 
coupled system of difference equations on the square lattice for w^^-' and w^'^\ we would get the same 
result. So the introduced potential equation p3p does not provide linearizable discrete equations. 

A discrete linearizable Burgers equation has been presented by Levi, Ragnisco and Bruschi considering 
Backlund transformation of the Burgers hierarchy UM, and by Heredero, Levi and Winternitz [QllTO] . In [M] 
we can find the discrete equation 

Um+l,n[p + '«m+l,™+l] - Um,n\p + Um+l,n] = 0, (61) 

and its Lax pair 

V'm,n+l — ^m,n^m,n^ ^'m+l.n — ' '4^ra,n- (^^) 

P + Um+l,n 



The linear PAE corresponding to ([CT]) is 



^m+l,n+l = "0™,™ -P'4^m+l,n- (63) 

Eqs. ([S^ suggest to rewrite ([55]) as 

^n+l,m/^n,m — ^m,^n(^^^,^7 ' * ' ^n,m+l/'^n,m — ^n^lni^^j^n: * ' * )' (^4) 

Eqs. p31 164p are transformable one into the other by defining = log(w„,m) and redefining appropri- 
ately the functions Sn.m and £n,m- However in doing so, if Wn,m satisfies a linear equation, this will not 



be the case for Vn.m- So the fact that the ansatz (|33p do not give rise to linearizable equations is not in 
contradiction with the fact that is linearizable. 
The compatibility of ([M)) implies 

If (|65p is constrained to be an equation on the square lattice, then we must have £n.m{un,m, ■ ■ ■) = 
£n}n{un,7n, Mn+i,m) and £n}n{un,m, ••■)== £n}n{un,m, Un,m+i)- Moreovcr with no loss of generality we can 

set ^n.m (Wn,m ; "^nj/n + l ) — '^n^m- 

Let us look for the symmetries of ([M]) . Applying the infinitesimal generator (|36p to the right hand 
equation in ([M)) we get 

■00,0 = V'o,o(wo,o), 00,0 = ■ (66) 

■^0,0 

Then the determining equation associated to the left hand equation in (j64p is given by 

V'i,o(^'i,o) = 77-^00.0 + -^-^01.0 "0,0 + '^^o^o'/'o,o(wo,o), (67) 

where the functions (pij are expressed in term of the functions "04,^ through (1661) . 

As we look for linearizable equations, from Theorem [5] it follows that we must have: 



V'o,o(wo,o) = XI ^o^o7o^d(wo,o), (68) 

where the discrete functions 'Wq \ satisfies a linear difference equation on the square. We can assume that 
the coefficient of \ is always different from zero so that we have 



In such a case the variables 



0) U) 

-^0,0' '"1.0 

„o-) 



of equations relating the functions 7o q(wo,o)5 J = 1; 2, with the function q 



and Wq\, j — 1,2, are independent and (|571) splits in three couples 

(1) 



7i,o'--o,o 



(2)c(l) _ 



7l, 0*^0,0 





(1) 



i9ui.o 



--0,0 



c-(l) 

-0,0 



du 
d£, 



(1) 

0,0 (1) 

— 7o>o,o 

0,0 

(1) 

0,0 (2) 

— 7o,>o,o 



(1) 

0,0 



duiA 



(1) 

"1,071,0 



(2) 

"i,o7i.o 



(1) 



9mo,o 



^(i)f^ (1) 
"^0,0 7o,i 



d£ 



9mi,i 
(1) 

0,0 



+ 


"--0,0 


+ 


"--0,0 



' (1) 

7o,o> 



' (2) 

7o,o> 



duin 



£> 



d£, 



(1) 

0,0 (2) 
— ^0,1 



d£ 



(1) . 



duo,o 



duin 



where vq.i = uo,oWo,o, ^^i,o = ^^0^0,0 and vi,i = ui,o£o^oWo,o due to ([Ml 



As £q^q is a function of uq.o and is a function of wo,o we get from (|70l [7T|) 



„o-) 



d£, 



(1) 
0,0 



9wi,o 

From dZl [721) we get 



'1,0 



7i,o 



'^o.o ~ '"o,oc-o,o 



(1)^ 

duo,o 



ai^h^ + ai-h?} ~ u,^,^^ di^h^ + di^h?} - "1,07^^ 



Omi,o 

while from ((74l [75)) we get 



^(1) 

'1,0 



'1,0 



c.(i)££ol 

"0,0 3 mo,o 



= '«2 



5ui 



When 1^ ^ 0, solving the equations for f (^'J from ([H [771 [75]) we get 

^(1) ^ ^2^^0,0 + ^1 ^ ^79) 
Ko - UlO 

With no loss of generality we can set 7''2) = and then e'--'-' = d'^' = /'^' = 0, j = 1, 2, w^^) — q and the 
equations ([751 [771 [751) are compatible if K2Koa'^^^ = ki6'^^^c*^^^. The resulting class of linearizable PAEs 
([7^ is an extension of the Burgers equation ([CT)) 



'1,0 



^(1) 



(2) 
^1,0 



(70) 
(71) 
(72) 
(73) 
(74) 
(75) 



(76) 



(77) 



(78) 



(kq - ui^o)(k2Uo,i + «i)uo,o - (ko - ui,i)(k2Wo,o + Ki)ui,o = 0, 



(80) 



which reduces to it when kq ^ 0, ki = 1 and K2 = 0. In ([501) all generality ki can be taken to be either 
or 1. Other two Burgers equations are obtained taking kq ^ 0, ki = and K2 ^ or kq = 0, ki = 1 
and K2 ^ 0. All the these three Burgers equations can be transformed to (1 + uo,o) ""1,0 = (1 + "0,1) ^0,0 
and we recover the results obtained in |24] . Moreover, if K2 7^ 0, ki = 1 and kq ^ 0, by the transformation 
"0.0 = Kq 1lZ°l g°'°+o!^, ' where {ej,Oj), j = 1,2 are arbitrary parameters, wq.o will satisfy the Hietarinta 
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equation 



^0,0 + 62 Mia + 02 



"1,0 + oi tto,i + ei ' 



(81) 




0, we must have fg q7}' 



= which has no nontrivial solution. 



5. Conclusions 



In this article we studied from the point of view of the symmetries the linearization through invert- 
ible or non-invertible transformations of nonlinear partial difference equations defined on a fixed non- 
transformable lattice. We find, in strict analogy to the continuous case, the linearizability conditions and 
apply them to the classification of linearizable PAEs on a lattice of four points in the plane. The results 
we obtain are compatible with the results obtained previously by requiring the existence of a linearizing 
transformation, however the classification turns out to be easier. 

Work is in progress to extend these results to the case of linearizable equations belonging to more 
general plane lattices and to nonlinear partial difference equations defined on transformable lattices. 
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In this Appendix we state and proof a theorem on the necessary and sufficient conditions for a partial 
difference equation to be linear. For the sake of simplicity of the presentation we limit our considerations 
to the case when the equation is defined on 4 points (see Fig.l), i.e. 
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Appendix 




(82) 



'n,m+l 



^n+l,m 



Figure 1 . Four points on a square 
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Theorem 6. Necessary and sufficient conditions for a discrete equation \8'^) to have a symmetry of 
infinitesimal generator X — (l)n,mdu„ ,„ , is that it is linear. The function <j)n.m must satisfy the following 
linear homogeneous discrete equation 

^(^4*n,m^ 4^n+l,mi 4*n,m+l^ 4^n+l,ni+lJ ~ 0n4-l,m+l ^n,ra^n,in ^n,m4^n,m+l ^n,m4^n+l,va — 0- (^3) 

Proof. It is almost immediate to prove that a linear partial difference equation defined on four lattice 
points (|82p has a symmetry (|27p where 4>n,m satisfies a homogeneous linear equation. To obtain this result 
it is just sufficient to solve the invariance condition 



prXf 



0. (84) 

£=0 



Not so easy is the proof that an equation which has such a symmetry (P7)) must be linear. In the full 
generality ([55)) can be rewritten as 

= F(n,m,(j) (85) 



and, by assumption, (1821) does not depend 

4^n,m and ()85p on Uji ^^i. 
Let us prolong the symmetry generator (|27p to all points contained in 

prX = (t>n,mdun^ + 4>n+l,mdun+i,m + m+l9it„,„+i + 4>n+l,m+ldun+i^^ + i- (86) 

The most generic equation (j82[) having the symmetry (j86p will be written in terms of its invariants 

£{n,m,K^,K2,Iu)=0, (87) 

with 

^1 = T- -r—^ -^2 = , A3 = ■ . (88) 

As ([57)) depends on n, m we can with no loss of generality replace the invariants ([55)) in ([57)) by the 
functions 

^ 4'n,m+l 7y 4'n+l,m /on\ 

-t^l — 'U'n,m+1 — Un,m—; , A 2 — Wn+l,m " U„,m — 7 , \o\) ) 



F(n,m,(j3 

K3 = Un+l,m+l — Un^m 7 



'-'n.ra 



Invariance of (1571) then requires = -g^frfi — = g^^^ = 0, i.e. 



d£ ( \ dE ( F F_^^ „ \ 

— ^ "n.m— 5 + 75 + — ^ 7 ^ "n.m = 0, (90) 

^ ( - J^) + ^ ( - = 0, (91) 

^ ( - + ^ ( - = 0. (92) 

As 7^ and ([5(11 IMl 1^ are a homogeneous system of algebraic equations, the determinant of the 

coefficients must be zero. Consequently the function F must satisfy the first order linear partial differential 
equation 

F — 0n,m^,0„,„ — '^n,m+l^,0„,„+i — 0n+l,m-F,0„+i,„ — (93) 

i.e. F is given by 

j-i , p I f. N f. 4'n+l,m <Pn,m+l fr,A\ 

F ^ (l)„^,nf{Lr), ^ = -7 ' ' 

where /(^, r) is an arbitrary function of its arguments. 
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Figure 2. Four points on a triangle. 



For F given by (|94l) the system ([TOl [Ml [M)) reduces to the foUowing two equations for 

oKj 

d£ d£ d£ d£ 

dKi ^OKs dK2 dKs 

whose solution is obtained by solving ([55]) on the characteristics 

£ = £{n,m,L), (96) 

r— _f — f ( „ '/'n.m + l \ _ r ( _ 0n+l,m \ 

V 0n,m / ■ V ■ 0n,m ^ 

Requiring that £ be independent of 0„,m, 4>n+i,m and 4>n./m+i we get /.^^ /.rr = /,Cr = i.e. 

^n+l,m+l ^n,m'^n,m ^n.m^n,m-\-l Cn.m^n+l.m- (^^) 

£ = is an (non autonomous, maybe transcendental) equation for L which, when solved, gives L = d„ 



where dn^m stands for the set of the zeros of the equation (in addition to n and m possibly dependent on 
a set of parameters). In conclusion Un.m must satisfy the linear equation 

^n+l,m+l ^n,m'^n,m ^n,m'^n,m+l ^7i.m^n+l.m ^7i,m ('^^) 

□ 



Remark 1. The proof of Theorem\^ does not depends on the position of the four lattice points considered 
in 12^. The same result is also valid if the four points are put on the triangle shown in Fig. 2, i.e. 

£ {n, m, U„_l,m, Un,m, W„+l,m , "n,m+l ) = 0. (100) 
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